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The very first data from MAX-lab in Lund, Sweden on pion photoproduction at 
threshold energies are presented. The decrease of the total ir + yield in 7 + 12 C, 
27 Al reactions below 200 MeV as well as the da/dQ cross-section data essentially 
follow the predictions of an intranuclear-cascade model with an attractive potential 
for the pion-nucleus interaction. However, d 2 a"/dOdT, cross-section data at 176 MeV 
show deviations which call for refinements of the model and possibly also for the 
inclusion of coherent pion-production mechanisms. 
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1 Introduction 



Studies of (7, it) reactions on nuclei provide information on the in-medium 
pion-nucleon (ttN) interaction [Tf2] . the properties of excited spin-isospin flip 
states of residual nuclei [314"] and the properties of few-body systems [5]. A 
proper theoretical description of photonuclear reactions and in particular of 
pion photoproduction on nuclei, is a prerequisite for calculations performed for 
the photodisintegration of ultra-high energy nuclei in the cosmic microwave 
background [6] and on nuclear fragmentation reactions induced by virtual 
photons [7]. 



Unfortunately, data on nuclear pion photoproduction are scarce in the near- 
threshold region because of lack of relevant photon beams. The recently up- 
graded nuclear-physics beamline at MAX-lab is one of only a handful of facil- 
ities worldwide which can now provide a photon beam of appropriate energy. 
This is done by colliding a pulse-stretched electron beam [8"f9"] with a thin 
radiator and momentum analyzing the post-bremsstrahlung electrons in one 
of two tagging spectrometers [10]. That said, for the early measurements re- 
ported on in this paper, the tagging spectrometers were not yet commissioned. 
Thus, in this untagged experiment, the photon-endpoint energy was 189 MeV 
and the average photon-beam energy was estimated to be 176 ± 2 MeV (see 
section 2.6). 



Some data on near-threshold pion production exist [3.4 ; JTT] but these come 
mainly from experiments with rather high pion detection thresholds. In this 
experiment, the two range telescopes allowed for a relatively low threshold 
for pion detection (8.4 and 8.7 MeV) which was achieved using a thin (3 
mm) first plastic detector and a compact telescope design. Unfortunately, the 
limited solid-angle coverage and the fact that the absolute normalization of 
the data reported on here utilizes (7,p) cross-sections from the literature result 
in large statistical and systematic uncertainties. 



The measured total cross-section data for tt + photoproduction on 2 H, 12 C 
and 27 Al, together with the corresponding differential (da/dQ) and double 
differential (d 2 a/dQdT) cross-section data have been compared with results 
from the RELDIS Monte Carlo model for photonuclear reactions p~2lTT3] . In 
this model, the quasi-deuteron absorption mechanism coexists with quasifree 
meson photoproduction on individual nucleons. 
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Fig. 1. Principle geometry and ir + detection process in the range telescopes (not to 
scale) . 



2 Experimental details 

2.1 Photon beam and targets 

The stretched electron beam from the MAX I ring [8] has long been used for 
photonuclear experiments at the Tagged-Photon Facility at MAX-lab. As pre- 
viously mentioned, in this, very first (7, tt) experiment, the tagging spectrome- 
ters [9] were not in operation. The electron beam was simply passing through a 
150 /im thick Al radiator which produced a photon beam with bremsstrahlung 
energy distribution in the interval - 189 MeV. These photons passed through 
a 19 mm diameter collimator before impinging on experimental targets of C, 
CD2 and Al with thickness 1-2 mm. All results on 7 + 2 H reactions have 
been determined from the data on C and CD 2 targets. 



2.2 Pion detectors 

The technique used to identify 7r + in this experiment was based on measuring 
the 26 ns 7r // decay in plastic scintillator range telescopes (Fig. 1). The effi- 
ciency for identifying pions depends here mainly - but not solely - on telescope 
geometry. It was important to minimize both the electronic noise through 
proper grounding and the room background through concrete shielding. After 
such precautions, we found that plastic scintillator telescopes of the CHIC 
design [H] were very effective at selecting pions. 
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Consequently, two telescopes of this type were designed for the experiment 
under the assumption of an effective upper limit in pion energy well below 50 
MeV. Each telescope was comprised of five NE102 plastic scintillators (Si - 
S 5 ), which were read out by Philips XP2020 photomultipliers (PMTs). Table 1 
presents the geometries chosen for the telescopes and the resulting energy bins 
for protons and pions that stop in each of the detector elements. Since pion 
identification required at least two energy signals (see below), the low-energy 
limit for pion detection was set essentially by the thickness of element S\. 
Including the minor energy loss that pions experience in target material, air 
and wrapping of the scintillators (which is different for the two telescopes), 
the detection thresholds for pions were 8.4 ± 0.1 MeV and 8.7 ± 0.1 MeV 
for the 30° and 90° telescopes respectively. The maximum pion energies each 
telescope could stop were 57.5 MeV and 49.9 MeV, respectively. The fact that 
one single pion was registered in each of the two last (S 5 ) detectors shows that 
the choice of telescope thicknesses was reasonable. 



Table 1 

Telescope geometry (th = thickness) together with the proton and pion energy bins 
(T p , T„-) for the telescopes. 





90° telescope 


30° telescope 


det. 


th 


area 


T P 


T. 


th 


area 


T P 


TV 




(mm) 


(mm 2 ) 


(MeV) 


(MeV) 


(mm) 


(mm 2 ) 


(MeV) 


(MeV) 


Si 


3 


60-60 


- 17.2 


- 8.7 


3 


60-60 


- 16.8 


- 8.4 


s 2 


5 


85-85 


17.2 - 29.1 


8.7 - 13.8 


5 


70-70 


16.8 - 28.9 


8.4 - 13.3 


s 3 


20 


90-90 


29.1 - 57.7 


13.8 - 26.2 


10 


80-80 


28.9 - 45.2 


13.3 - 20.4 


S4 


30 


95-95 


57.7 - 86.4 


26.2 - 38.8 


20 


90-90 


45.2 - 68.3 


20.4 - 30.2 


s 5 


50 


100-100 


86.4 - 122.6 


38.8 - 57.5 


50 


100-100 


68.3 - 109.2 


30.2 - 49.9 



2.3 Electronics and data acquisition 



In order to identify 7r + , the PMT anode signal of each detector was split and 
sent to two independent charge integrating ADCs with individually adjustable 
gates. The "prompt-gated" ADC had a 100 ns long gate which was open for 
integration 5 ns before the analogue PMT signal appeared. The "delay-gated" 
ADC, also had a 100 ns long gate, which was opened when the largest ampli- 
tude of the analogue pulse was reached. This opening time was individually 
adjusted from the proton pulses, easy to recognize with beam on target. 
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For each master trigger which required signals in the first two detector ele- 
ments, S1*S2, the hit-pattern was registered and stored in ROOT-tree files 
together with all the ADC and TDC information. A local, VME based data 
acquisition system was used for this storage. 

2.4 Pion identification 

Figure 2 presents typical on-line data. The upper plot shows the prompt AE 
- E correlation (in this case, S2 vs S3 in the 30° telescope) after the ADC 
and TDC patterns have been used to identify the detector element in which 
the pion stopped. No other constraints have been imposed. The existence of 
three groups of particles - protons, pions and e-7 background - is obvious. It 
is important to note that the PMT gains were set so that all protons were 
registered (see section 2.6). 

The pion band contains both tt + and 7r~. 7r + were then selected in the E prompt 
vs. Edelay plots for the "stop" detector (Fig. 2 - mid). The muon from the 
7r + — >• fi + u decay (nearly 100% branching ratio) can be detected with a certain 
efficiency (see section 2.5). The 7r~ is normally absorbed by a carbon atom of 
the scintillator and it thus appears in the E prompt - Edeiay plot either in the 
lower band if only neutral particles are emitted from the disintegrating carbon 
nucleus or in a random background position if charged particles are emitted. 
Selecting the upper band in Fig. 2 (mid) thus provides clean tt + identification 
as evidenced by Fig. 2 (lower), which shows only those events originally present 
in Fig. 2 (upper) which pass the n + identification. 

Pion-energy intervals were defined from the energy-range intervals that stopped 
pions have in the detectors (Table 2). Polynomial fits, to the energy distribu- 
tions, using the detector half thicknesses were then used to determine the final 
energies. This resulted in systematic uncertainties of ~ ±1 MeV in the pion 
energies. Note that linear scintillator response functions have been used when 
pion detection thresholds are artificially introduced in software. 

2. 5 Determination of the ir + detection efficiency 

The number of n + obtained in the three-step analysis process described above 
was corrected for the following processes: 

i) pion decay in flight. 

ii) nuclear reactions in the target and detector material. 

iii) sliding trajectories (detector geometry). 

iv) ji + detection efficiency. 
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Fig. 2. (Upper) Prompt AE - E plot for particles stopping in the S3 detector of the 
90° telescope. (Mid) Delayed vs. prompt ADC signals for events lying in the pion 
band in the upper plot. (Lower) Prompt AE - E plot revisited, now for those events 
defined as ir + in the middle plot. 



The first three processes make pions undetectable, while correction iv) ac- 
counts for losses of n + in the E prompt - E de i ay plot when muons leave the stop 
detector before decaying (/i + — ► e + w). In the efficiency calculations, it was 
assumed that half of the muon energy (0.5 ■ 4.2 MeV) must be deposited in 
the detector for the pion to appear in the proper band in the E prompt - E de i ay 
plot. All effects have been calculated analytically following the procedure in 
Ref. [H] and also by Monte-Carlo simulations using GEANT4 [15J. Table 2 
shows the results. Note that the analytically determined, total corrections 
(ranging from 24 to 44%) are systematically lower than the results from the 
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Monte-Carlo simulations (32 - 69%). This difference is interpreted to result 
from approximations in the analytic scattering kinematics calculations. The 
larger correction factors resulting from the simulations were thus used to ob- 
tain the true number of photopions. The differences in the corrections were the 
main contributions to the systematic uncertainties in the detection efficiency 
(see section 2.7). 

Table 2 

Total correction factors for detector efficiency in the elements S2 to S5 calculated 
analytically and simulated with GEANT4 (see text for details). 





a (90°) 


telescope 


P (30°) telescope 


Det. 


Analytic corr. M-C corr. 


Analytic corr. 


M-C corr. 


s 2 


1.44 


1.69 


1.44 


1.69 


S3 


1.31 


1.48 


1.33 


1.54 


S4 


1.28 


1.39 


1.30 


1.42 


s 5 


1.24 


1.32 


1.26 


1.37 



2.6 Normalization and average beam energy 



Because of the low data acquisition rate, <C T10 3 events/s, no proton rejection 
was necessary. Consequently, the registered protons could be used for absolute 
normalization. A mixture of inclusive and semi-exclusive (j,p) cross-sections 
have been utilized for this normalization [T61I17II18II19II20I. 



The fundamental, triple differential cross-section that was extracted from this 
experiment is, 



-7T,p 



where T^p and Q WjP are the kinetic energy and emission angle of protons and 
pions respectively. C^p are the normalization constants (which depend on tar- 
get mass number, thickness and density) and $ (E 7 ) is the bremsstrahlung 
spectrum of incident photon energies. N^p is the efficiency-corrected pion(proton) 
yield. AQ is the solid angle subtended by the detector. The most extensive 
(7,p) data were found for photon collisions on 12 C. Total, inclusive cross- 
sections could in this case be found at energies > 200 MeV [16] while total 
cross-sections at lower energies must be estimated from missing energy biased 
data [17J. The normalization constant C p was obtained from the number of 
protons registered in all of the stop detectors, J2^ 9 =o Y^'=i7 N P {E 1: 6 P , T p ) 
and the empirical triple differential cross-section integrated over the proper 
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intervals in E 7 and T p . In practice, a polynomial fit to the function f(E 7 ) = 
/o 189 $ (£ 7 ) ■ jg" dEj fn pdTp ( E ^ Q P ,T p )dE^dT p was extracted. 

If the pion triple differential cross-section had been known, C n = k ■ C p 
could be determined in a similar manner (and ideally k = 1). Instead, an 
ansatz was made by extrapolating a polynomial fit to the Fissum et al. data 
[IT] (photon energy region 184 - 213 MeV) to an energy of 160 MeV, be- 
low which the pion contribution was neglected. was then calculated in 
the same way as C p . Thus, a polynomial function was fitted to g(E 7 ) = 
Jo 189 $ (£ 7 ) • /g* — dEi f n a wdTn (£ 7 , B^, TjdE^dT*. The number of pions in- 
cludes here only those with T n > 17 MeV, efficiency corrected accordingly 
(Table 1). As mentioned before (section 2.4), linear scintillator response func- 
tions were assumed when fractions of "stop" detector energy bins were intro- 
duced. In a typical example of this procedure the determination of absolute 
differential cross-sections for 7 + 12 C at 90° k was determined to 0.79 with an 
uncertainty of 15% (see section 2.7). 

The same procedure was carried out for the 7 + 27 A1 and 7 + 2 H reactions. 
Because of the limited number of " normalizing" data, available for 7 + 27 Al, 
also 7 + 40 Ar data were used to guide the A T dependence of C n . 

The ansatz for the g(E 7 ) function resulted in < E 1 > = 178±2 MeV for pion 
producing events. Introducing the data from this work at 178 MeV photon 
energy, resulted in a decrease of < E 1 > to 176±2 MeV. Further recursive 
iterations along this line changed the value by < 1 MeV. The < E 1 > was 
consequently fixed at 176±2 MeV. Within the limits of uncertainty, the aver- 
age photon energy was the same for all three 7 + 27 Al, 12 C, 2 H reactions. 



2. 7 Systematic and statistical uncertainties 

Statistical uncertainties only are presented with the data in Figs. 3-7 with two 
exceptions - the two upper points (triangles) in Fig. 4. These contain the yield 
of pions with energy below the detector cutoff, which has been extrapolated 
both from data and RELDIS calculations. The extrapolated yields contribute 
with 30% and 50% of the differential cross-sections at 30° and 90° lab angle 
respectively. The uncertainties of these yields, which are introduced in the 
error bars of the two points in Fig. 4, are set to 30 %, based on the difference 
between the empirical and RELDIS extrapolations. 

Three major sources of systematic uncertainties in eq. (1) were identified. 
First, the systematic uncertainty in the detector efficiency has been taken 
as the difference in the results obtained from the analytic and Monte-Carlo 
calculations. This gives a pion-energy dependent uncertainty ranging from 
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7% at high energy to 17% at low energy (Table 2). The second source of 
uncertainty comes from the manner in which the cuts identifying the pion yield 
was applied. This uncertainty was estimated to be ~ 9%. The third systematic 
uncertainty enters via C p through its dependence on the (j,p) cross-section 
data taken from the literature. By comparing the pion yields from the 7 + 
12 C reaction that are obtained with and without the use of the data from Ref. 
[T6] to estimate C p , this uncertainty was determined to be 12%. Other sources 
of uncertainty were negligible. This resulted in total systematic uncertainties 
for the 7 + 12 C reaction that range from 17% at the highest pion energy to 
23% at the lowest. The uncertainties are somewhat larger (up to 28%) for the 
other two reactions. 



3 The RELDIS model for photonuclear reactions 



Below pion production threshold, at E 7 ~ 140 MeV, the de Broglie wave- 
length A is comparable to the distance between nucleons in nuclei and photon 
absorption by a quasi-deuteron is the main reaction mechanism. At higher 
energies (as A becomes comparable to the nucleon radius), photons interact 
mainly with single nucleons, thereby exciting baryon resonances and produc- 
ing mesons. The RELDIS model takes into account these two competing chan- 
nels. Here, the input to the code relevant for photon energies close to the pion 
production threshold is described. More details can be found in Refs. [T2"][T3"] . 



3.1 Calculation of quasi-deuteron absorption 



The two-nucleon photoabsorption cross-section on a heavy nucleus is 
taken from the quasi-deuteron model of Levinger |21j . as modified in Ref. |22j, 

a*]? = kZ(l - ZjA)aT h . (2) 



Here a e d is the meson-exchange part of the cross-section ad for deuteron 
photodisintegration, 7c/ — > np [23], A and Z are the mass and charge numbers 
of the target nucleus and k « 11 is an empirical constant taken from the 
analysis of Ref. [22J . 

The cross-section er^f decreases strongly with photon energy. Nevertheless, 
the two-nucleon absorption mechanism competes noticeably with the single- 
nucleon absorption up to E y ~ 0.5 GeV. 
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3.2 Simulation of pion photoproduction on nucleons in a nucleus 

The main single-nucleon photoabsorption mechanism is meson production. 
The two-body channel — ► ttN dominates up to E 7 ~ 0.5 GeV. This 
process was calculated in the framework of the phenomenological approach 
of Ref. [2"H25] . In this approach, the pion photoproduction amplitude con- 
tains Breit-Wigner resonant terms, Born terms and a weakly energy- dependent 
"background" contribution. The latter was used as an adjustable parameter. 
Masses and widths of the resonances were taken from 7rN-scattering data, 
with their amplitudes taken as free parameters. The excitation of six different 
baryon resonances were considered with the A(1232), iV*(1520) and iV*(1680) 
resonances the most important of them. 

Tables of total reaction cross-section data together with photopion angular 
distributions calculated according to Refs. [24.25J were taken from Ref. [26J. 

Both the total and the differential cross-sections of the 7p — > 7i + n, 772 — > n~p 
and 7p — ► ir°p processes are well described by the jN Monte Carlo event 
generator [12] used by the RELDIS model. Presently, the jN event generator 
used in the RELDIS model extends the event generator by Corvisiero et al. [26] 
to higher photon energies (up to 10 GeV) and multiple pion production, as 
demonstrated in Refs. [T2"fT3] . 

3.3 Secondary interactions of photohadrons 

Hadrons produced in a primary 'yN or 7c? — > np interaction initiate a cascade 
of successive hadron-nucleon collisions inside the target nucleus during the 
cascade stage of the photonuclear reaction. Calculations are based on a Monte 
Carlo technique used to solve the equation that describes hadron transport 
in the nuclear medium. The target nucleus is considered to be a mixture of 
degenerate Fermi gases of neutrons and protons in a spherical potential well 
with a diffuse boundary. By using the effective real potentials for nucleons 
and pions, the influence of the nuclear medium on cascade particles is taken 
into account. It should be stressed that this potential is taken to be constant 
although it should depend on the nuclear density as well as on the pion kinetic 
energy. However, it is beyond the scope of the present experiment to provide 
the very precise data necessary to tune these dependences. 

The momentum distributions of nucleons in the nuclei are calculated in the 
local-density approximation of the Fermi-gas model. The distribution of nu- 
clear density is approximated by a set of step functions for the nuclear radius. 
The Coulomb potentials for charged cascade particles are calculated for each 
density zone. 
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The cross-sections for pion interactions such as ttN — > ttN, tt(NN) — > NN, 
ttN — > ttttN etc, as well as nucleon-induced processes, A^ A" — > iViV, A^ A" — > 
ttNN, ,. . . etc in the nuclear medium are assumed to be the same as in vacuum, 
except that the Pauli principle prohibits the transition of the cascade nucleons 
into states already occupied. 



4 Experimental data and comparison to RELDIS calculations 

Data on 7r + cross-sections were measured in the TV interval(s) 8.7(8.4) - 
57.5(49.9) MeV. No pions with higher energies are expected (section 2.2). 
For simplicity, these pions are subsequently denoted as T^ > 9 MeV pions. 
Due to limited statistics, data points are presented one per "stop" detector 
(Table 1), all at beam energy E 7 = 176 (± 2) MeV (recall section 2.6). 

In the near-threshold RELDIS simulations used for comparison to the data 
in Figs. 3-7, the only process for 7r + production is 7p — > 7T + n. Photopions 
can then face elastic scattering from nucleons in the target nucleus, charge- 
exchange reactions (such as 7f + n — > ir°p) or absorption on two nucleons. The 
key parameter in these simulations is the real part of the pion potential V v 
which is assumed to be attractive, < 0. This potential is taken as an 
empirical parameter, independent of the pion energy. Its value has been derived 
from pion production data reasonably close to threshold. The imaginary part 
of the pion potential is simulated by taking into account the pion absorption 
reaction tc(NN) -> NN. 

Because of the limited statistics and large systematic uncertainties, it is dif- 
ficult to compare the double differential cross-section data to the predictions 
from RELDIS, which themselves are very sensitive to small changes in the 
input parameters. Therefore, total (energy and angle integrated) cross-section 
data are presented and discussed first. Least squares fits to Boltzmann func- 
tions (energy) and polynomial functions (angular) have been introduced to 
carry out these integrations. 

4-1 Total yield o/7r + 

Data from Ref. [TT] demonstrate a falloff of a n + with decreasing E 7 towards 
the effective reaction threshold. The measured total cross-section data for 
pions with T w > 17 MeV (the cutoff in Ref. [UJ) from this experiment and 
from selected parts of Ref. [11] are shown in Table 3. The average systematic 
uncerainty in the 1 H and 12 C data from Ref. [TIj was ~ 17%, i.e. somewhat 
less than in the present experiment (see Table 3). 
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Table 3 

Total cross-section data for ir + photoproduction with T^ > 17 MeV from this work 
and from Ref. [11] . Statistical uncertainties are shown. 



a(pb) 


E 7 (MeV) 


X H 2 H 


12 C 


27 Al 40 Ca a 


176 (this work) 


0.8±0.2 b 


7.2±1.2 c 


43±7 d 


184 p] 




31.1±2.0 


131.8±7.8 


194 p] 


82.0±3.7 


64.4±2.8 


177.8±11.7 


204 pi] 


92.3±4.2 


107.8±3.8 


354.3±13.0 


213 [IE] 


103.7±4.7 


160.7±4.5 


498.1±15.0 



a Systematic error associated with the background 8.8 jjh [TT] 



b Systematic error 0.2 /A>, c Systematic error 1.4 fib, d Systematic error 12 /A> 

The process that RELDIS uses to produce pions should result in a n + scaling 
with Z target . Thus, <7 W + I Ztarget vs. E 7 is plotted in Fig 3. The Fissum et al. 
data show that heavy nuclei are less efficient in photoproducing pions. This 
may be due to the enhanced reabsorption of tt + . The strongest deviation from 
^target scaling appears actually in the 7 + 1 H cross-section, which is 5 - 10 
times more efficient in producing ir + than the 7 + nucleus reactions [TT] . 

The new data shown in Fig. 3 exhibit an enhanced falloff of pion production 
in the 7 + 12 C reaction, consistent with the predictions of the RELDIS model. 
There is of course a shift from the (7, N) channel to the quasi-deuteron ab- 
sorption channel in this energy region but on the other hand, the 7 + 27 A1 data 
do not show this enhanced falloff. This is in contradiction to the RELDIS re- 
sults and the deviation is well outside statistical and systematic uncertainties. 
RELDIS results also exhibit a steeper falloff for heavier nuclei (note that Fig. 
3 shows 7 + 40 Ca calculations). The enhanced production in heavier nuclei 
exhibited by the data taken in this experiment, may be due to an underesti- 
mation of the tail of the internal momentum distribution of the nucleons or a 
signal that coherent processes should be included in the model. 

Another remark to Fig. 3 is that the ir + cross-section in the 7 + 2 H reaction 
does not follow the elementary 7 + p trend, but rather the trend established 
by the 7 + nucleus reactions. It thus appears as if the presence of one single 
pn pair leads to the normal competition between the single nucleon and quasi- 
deuteron processes. 

It should finally be stressed that in this near-threshold energy domain, the 
available phase-space for pion production becomes increasingly limited and a 
high detection threshold (e.g. 17 MeV at 90°) makes this restriction even more 
important. 
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Fig. 3. Total yield of tt + with TV > 17 MeV (the pion-detection threshold from Ref. 
|llj ) as a function of E~. Points represent data (this work and Ref. [11]) according 
to the legend and curves represent RELDIS calculations with a real pion potential 
of -20 MeV for the 7 + 12 C (solid) and 7 + 40 Ca (dashed) reactions. Systematic 
uncertainties of the measurements are given in Table 3. See text for details. 



4-2 Angular distributions of produced it 



In Fig. 3, the falloff of ir + photoproduction at energies < 180 MeV was noted. 
Fig. 4 shows that the differential cross-section, da/dQ, decreases by a factor 
4 - 7 if the detector limit is raised from 9 MeV (this experiment) to 17 MeV 
(detection threshold in Ref. [II]). The extrapolated increase in pion yield if 
the detection threshold could be lowered to MeV, is a factor of 1.5 - 2. The 
trend exhibited by the data is reasonably well reproduced by RELDIS, with 
V,r = -20 MeV. Systematic uncertainties associated with the data points are 
are between 17% and 23% and the angular acceptance is ± 5.4° and ± 7.3° 
for the forward and 90° telescopes respectively. 

This is also essentially true for the 7 + 27 Al reaction, as shown in Fig. 5. 
The excess pion yield was here essentially found in the forward hemisphere. 
This fact goes hand-in-hand with the assumption of an extended internal 
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Fig. 4. Angular distributions of ir + with TV > 17 MeV (open points), > 9 MeV 
(squares) and > MeV (triangles), emitted in 7 + 12 C reactions at 176 MeV. 
Systematic uncertainties of the measurements are between 17 and 23%. Curves 
represent RELDIS calculations with a real pion potential of -20 MeV. 



momentum distribution since the coupling of the relative momentum vector 
to the internal momentum vector produces a forward-focusing effect. However, 
there are other tentative explanations such as too strong reabsorption in the 
model that are plausible. Finally, we note again that the 7 + 2 H reaction 
behaves more like the reaction on heavier nuclei than the elementary process. 



4-3 Double- differential cross-section 



The conclusion from the da/dfl data of tt + in 7 + 12 C reactions (Fig. 4) was 
a slight underprediction by RELDIS with V,,- = -20 MeV which, however, 
rather tends to become an overprediction when extrapolating the data to TV 
= MeV. The behaviour of the double differential cross-section as a func- 
tion of pion kinetic energy (Fig. 6) at 30° and 90° makes the overestimation 
even more plausible. The obvious overprediction for low TV is to some extent 
compensated by more extended high energy tails in the data. 
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©,ab ( de 9) 

Fig. 5. Angular distributions of tt + with T n > 9 MeV in 7 + A (see legend) reactions 
at 176 MeV. The systematic uncertainties are between 17% and 28%. The angular 
acceptance is ± 5.4° and ± 7.3° for the 30° and 90° respectively. Curves represent 
RELDIS calculations for A = 12 C (dashed) and A = 27 Al (solid) with a real pion 
potential of -20 MeV. 

The same tendencies are more pronounced for the 7 + 27 A1 data (Fig. 7). In 
this case, RELDIS calculations exhibit no high-energy tail and the expected 
maximum in a Coulomb shifted spectrum of ir + from the one-nucleon scatter- 
ing process seems to fall at lower energies in the calculations. All calculated 
curves in Fig. 6 and 7 have a steep rise in the second bin, corresponding to 
an effective Coulomb barrier of ~1 MeV. This appears to be too low, at least 
for the 27 Al case, where a more realistic Coulomb shift of 3 - 4 MeV would 
improve the comparison. The choice of the optimal pion potential may suffer 
from this mismatch. More data, especially on reactions with heavy targets, 
are however needed to set this question. At present it appears reasonable to 
believe that the potential falls in the interval -20 MeV to -15 MeV for light 
nuclei, with the lower value more representing 27 A1 and the higher (-15 MeV) 
representing 12 C. 

The form of the 30° energy spectrum of tt + from 7 + 2 H reactions (Fig. 7) 
is somewhat different in that the low energy point (~11 MeV) has a low 
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T (MeV) 

Fig. 6. Double-differential cross-section data for tt + photoproduction with T n > 9 
MeV in 7 + 12 C reactions at 176 MeV. The detector angle is 30° (upper panel) 
and 90° (lower panel). Note the experimental points at T n > 20 MeV. Systematic 
uncertainties are between 20% and 25% and the uncertainty of the energy positions 
typically 1 MeV. Curves represent RELDIS calculations with a real pion potential 
of -20 MeV (dashed) or -15 MeV (solid). 
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Fig. 7. Double-differential cross-section data for ir + photoproduction detected at 
30°, with TV > 9 MeV in 7 + 2 H, 12 C and 27 Al (see legend) reactions at 176 MeV. 
Systematic uncertainties are between 22% and 28% and the uncertainy in energy 
position typically 1 MeV. Curves represent RELDIS calculations for Al and C re- 
actions with a real pion potential of -20 MeV. 



d 2 o~/dfidT value. Even if this fact is significant, it must be stressed that the 
H 2 data are derived as a difference between the CH 2 target and C target data, 
which makes the systematic uncertainties larger, 26 - 28%, than for the other 
targets. 



5 Conclusions 



Near-threshold (7,71*) reactions in light nuclei may be described by an intranu- 
clear cascade model like RELDIS. This means that the basic processes in 
RELDIS at these energies namely quasi-deuteron absorption and photon ab- 
sorption on a single bound nucleon, should be dominant. Chosing depth of 
the real pion potential of -20 to -15 MeV results in RELDIS predictions which 
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match the data. Our approach to pion-nucleus reactions is similar to the one 
adopted in the Liege intranuclear cascade model [27J, which was recently re- 
considered in theoretical investigations of the average potential energy felt by 
a pion inside the nucleus. Our finding of V w = -20 MeV agrees well with the 
values found in the literature [27], with a surface component of V n between 
-5 and -25 MeV. Since low-energy pions, produced in the central part of a 
nucleus, are more easily absorbed than those produced close to the nuclear 
surface, it is relevant to use this surface component of V n for the comparison. 

The shapes of the experimental 7r + angular and (in particular) energy distri- 
butions differ from RELDIS predictions. Furthermore, the total tc + yield in 
heavier nuclei is underestimated by the code by a large factor. Model particu- 
lars, such as the effective Coulomb potential, the exact form of the degenerate 
nucleon momentum distribution and the energy and density dependence of 
the pion potential could certainly be tuned to improve the comparison. The 
discrepancies also hint that multinucleon channels can play a role not only in 
pion absorption but also in pion production processes 

7r + cross-sections per proton are much higher in the elementary 7 + p reac- 
tion than in 7 + nucleus reactions close to threshold. The fact that the 7 
+ 2 H data follow the 7 + nucleus trend indicates that the balance between 
quasi-deuteron absorption and single-nucleon absorption is achieved already 
for deuteron reactions. In fact, all data from this experiment indicate that the 
effects of the nuclear environment are small for nuclei with A < 30. 

We look forward to studying the details of the pion production mechanism 
in forthcoming tagged (j,ir) experiments at MAX-lab at energies up to ~220 
MeV in the very near future. 
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